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ABSTRACT 



Aims. We have investigated the nature of the variability of CHS 7797, an unusual periodic variable in the Orion Nebula Cluster. 
Methods. An extensive I-band photometric data set of CHS 7797 was compiled between 2004-2010 using various telescopes. Further 
optical data have been collected in R and z' bands. In addition, simultaneous observations of the ONC region including CHS 7797 
were performed in the I, J, K s & IRAC 3.6 and 4.5 [im bands over a time interval of as 40 d. 

Results. CHS 7797 shows an unusual large-amplitude variation of a 1.7 mag in the R, I, and z' bands with a period 17.786 ± 0.03 d 
(FAP = 1jc10 -15 %). The amplitude of the brightness modulation decreases only slightly at longer wavelengths. The star is faint 
during « 2/3 of the period and the shape of the phased light-curves for the seven different observing seasons shows minor changes 
and small-amplitude variations. Interestingly, there are no significant colour-flux correlations for A < 2 /mi, while the object becomes 
redder when fainter at longer wavelengths. CHS 7797 has a spectral type of M6 and an estimated mass between 0.04-0.1 M Q . 
Conclusions. The analysis of the data suggests that the periodic variability of CHS 7797 is most probably caused by an orbital motion. 
Variability as a result of rotational brightness modulation by a hot spot is excluded by the lack of any colour-brightness correlation 
in the optical. The latter indicates that CHS 7797 is most probably occulted by circumstellar matter in which grains have grown from 
typical 0.1 fim to ~ 1-2 yum sizes. We discuss two possible scenarios in which CHS 7797 is periodically eclipsed by structures in a 
disc, namely that CHS 7797 is a single object with a circumstellar disc, or that CHS 7797 is a binary system, similar to KH 15D, in 
which an inclined circumbinary disc is responsible of the variability. Possible reasons for the typical 0.3 mag variations in I-band at a 
given phase are discussed. 

Key words. Stars: low-mass, pre-main sequence - Stars: rotation, starspots, binary - (Stars:) circumstellar matter - Technique: photo- 
metric 



1. Introduction 

The photometric variability of T Tauri stars (TTS) is attributed 
in many cases to magnetically induced cool starspots and/or 
magnetically channelled variable accretion flows, which pro- 
duce hot spots at the base of the magnetic channels. During 
the last 10-15 years, various extensive CCD monitoring pro- 
grams of young clusters like NGC 2264 and the Orion Nebular 
Cluster (ONC) have been carried out. These monitoring pro- 
grams have significantly increased our knowledge of TTS vari- 
ability and rotation, and have revealed thousands of mostly low- 
amplitude periodic pre-main sequence (PMS) varia bles and a 
similar number of irregular variable PMS stars (e.g. | Herbst et all 
l2002t lLamm et all l2004t lRodrfguez-Ledesmaetal.ll2009l) . As 
reviewed by iHerbst et alj d2007l) . and based on many of these 
studies, it is possible to distinguish five types of PMS variabil- 



ity. The different types have known properties such as typical 
amplitude range of the variability at different wavelengths. For 
example, amplitudes at optical wavelengths of periodic bright- 
ness modulation, due to hot spots and/or magnetically chan- 
nelled accretion, can be a factor 2 to 5 larger than the peri- 
odic variations believed to be cau sed by cool spots, w hich can- 
not be larger than 0.5 mag in I (IHerbst et all 11991 . Another 
relevant difference between these two cases is that cool spots 
have been observed to be stable over thousands of cycles, while 
accretion-related variability is much more unstable and irregular, 
with the observed periodicity often only surviving several tens 
of cycles. It is important to note that rotational properties and 
accretion-relate d variability have been observed also in substel- 
lar objects (e.g.lScholz & Eisl6ffeil l2004al: iMohantv et~al]|2005l: 
Rodrigue z-Ledesma et al.ll2009l) . 
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Fig. 1. 1-band images of CHS 7797 during maximum and minimum brightness obtained with WFI at the 2.2 m telescope on La Silla, 
Chile. North is up and east is left. The region shown is 2 'x2 ' in size. 



Because most of these variability studies involve a large 
number of objects, they are not only valuable for characteris- 
ing variability, but also for detecting interesting and rare objects. 
One striking example is KH 15D, a unique 48 d periodic vari- 
able in NGC 2264, with extremely de ep (^3.5 mag) minima. I t 
was first recognized a s interesting by Kea rns & Herbstl dl998l) . 
lHamilton etall (l200ll) reported that KH 15D is a K7 PMS mem- 
ber of NG C 2264 (^ 2-4 Myr) with a mass of 0.5-1 M Q , and to- 
gether with lHerbst et al.l (120021) proposed that the star was being 
eclipsed by circumstellar material. Further intense investigations 
and, in particular, radial velocity studies led to a binary model 
for KH 15D in which the system is surrounded by a nearly edge- 
on and slowly precessing ci rcumbinary disc, to which the binary 
system is slightly inclined dWinn et alJl2006t iHerbst et al.l f2008. 
and references therein). Over the past 15 years the circumbinary 
disc has completely occulted the orbit of component B and has 
allowed us to see only a diminishing fraction of the orbit of star 
A. It is the appearing and disappearing of star A behind the disc 
that causes the strong light modulation. 

Other TTS binaries are expected to have similar character- 
istics, but KH15D is not only special because of its particular 
geometric orientation and the precessing circumbinary disc but 
also beca use of the indicatio ns of significant grain growth within 
the disc dHerbst et al.ll2008l) . This binary is also interesting be- 
cause it is a jet-driving source, with the jet being most probably a 
"common product" of the wh ole binary system or the circumbi- 
nary disc (Mu ndt et aDl20 10). Although thousands of PMS stars 
have been photometrically monitored in the past years, only few 
objects show a variable behav iour that resembles KH 15D , i.e. 
YLW 16A and WL4 in o Qph dPlaychan et al.l l2010. 20081, and 
V718Per in IC 348 dGrinin et al.ll2008l) . In the case of WL4 
and YLW 16 A, the authors proposed the existence of a new 
class of periodic disc-eclipsing binaries to explain the KH 15D 
type variability, in which a third body in these systems is re- 
sponsible for the warping of th e inner circumbinary disc needed 
to produce the eclipses, while iGrinin et al.l d2008l) argued that 
V718 Per is most probably a single star surrounded by an edge- 
on circumstellar disc with an irregular mass distribution at the 
inner edge of the disc, which causes the observed periodic vari- 
ability. AATau-like objects show also high amplitude variabil- 
ity, although this variability is found to be only quasi-periodic 



(lAlencar et al .1 l2010h . As described in lAlencar eta I] d2010h . in 
the fairly common AATau-like objects the quasi-periodic vari- 
abilit y is attributed to a magnetically controlled inner disc warp. 

Rodriguez-Ledesm a et al.l d2009) carried out a photometric 
monitoring in the ONC during 2004 with the primary goal of 
deriving rotational periods for a large sample of very low-mass 
PMS stars and brown dwarfs. Several of the objects in this sam- 
ple show unusual li ght curves: CHS 7797 (i.e. star No. 7797 in 
Carp enter et al.l (1200 ll) ) stuck out particularly because of its un- 
usual large-amplitude modulation of « 1.7 mag in I and its peri- 
odic variability. Hoping that this newly discovered periodic vari- 
able might have a similar potential to KH 15D, which has pro- 
vided considerable insight into the physics of PMS binary sys- 
tems and their circumbinary discs, we have carried out a photo- 
metric follow-up campaign, including multicolour data collected 
with various telescopes/instruments. 

This paper is organised as follows: in Section 2, we present 
the optical and near-infrared (NIR) data sets, while in Section 3 
we describe the photometric data evaluation. Section 4 describes 
the time series analysis procedure. We analyse colour changes 
in Section lBTTI and estimate luminosity and mass ranges for 
CHS 7797 in Section l5T2l A detailed discussion of our results 
is given in Section|6] while final conclusions are presented in 
Section 7. 

2. Optical and near-infrared data set 

CHS 7797 (7M(2000) = 05 : 35 : 08.7 D£C(2000) = 
-05 : 31 : 27.3) is located *9' south-west of the centre of 
the ONC . Fig.Q] shows I-band ima ges of CHS 7797 taken at 
maxim um and minimum brightness. Rodriguez -Ledesma et al.l 
(2009) were only able to derive an approximate period of 
18-20d from their data. This initial discovery triggered our 
photometric follow-up observations at various telescopes dur- 
ing the years 2009 and 2010. Observations were carried out 
from November 2009 to November 2010 with the 2.2m tele- 
scope on LaSilla/Chile, the 1.2m MONET/N Telescope at 
McDonald Observatory/USA, and with the 1.5m telescope at 
Mount Maidanak Observatory (MMO) in Uzbekistan. In addi- 
tion, we have also made use of the extensive archival data of the 
ONC (including the region around CHS 7797) collected between 
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Fig. 2. 1-band light-curves from the 2004 to the 2010 observing 
season. Observations from the different sites are shown with dif- 
ferent symbols (see Table 1 and text for details). 

Table 1. Overview on the obtained photometric observations. 



Site 


Observing^, 


Bands 


Number of 




seasons 




observations 


La Silla 2.2m (WFI) 


A,F 


I 


106 


Maidanak 1.5m 


F,G 


LR,z' 


233 


HCT and VBT (2 and 2.3 m) 


A-F 


LR 


346 


MONET 1.2m 


F 


I 


24 


UKIRT (WFCAM) 


F 


J 


59 


CFHT (WIRCam) 


F 


K s 


10 


Spitzer (IRAC) 


F 


3.6pm, 4.5/rni 


88 



" Observing seasons are denned as follows: 
A=2004, B=2005, C=2006/2007, D=2007/2008, E=2008/2009, 
F=2009/2010, G=2010 



December 4 2009, which were taken at this overlapping time in- 
terval by chance (see Fig. [5). This rich multi-wavelength data 
set allows us to put constraints on the physical characteristics 
of CHS 7797. Particularly for young and highly variable objects, 
such simultaneous data are of great importance for understand- 
ing their nature. 



3. Photometry 

For the observations taken with the 2.2m, MMO, and 
MONET/N, aperture photometry was performed using standard 
IRAF tasks. A set of known non-variable or very low vari- 
able stars close to CHS 7797 were selected to perform relative 
photometry between these secondary standards and CHS 7797. 
These non-variable stars were selected from the detai led moni- 
toring of the field by Rodrig uez-Ledesma et al.l J2009). With the 
MMO observations of Landolt standard star fields and SDSS im- 
ages the absolute magnitudes of these secondary standards in R, 
I and z' could be determined. A table with these magnitudes is 
given in the appendix. For the observat i ons ta ken at HCT and 
VBT we refer the reader to Parihar et al. 

Fig.|2]shows the I-band light-curve of the whole data set cov- 
ering the seven observing seasons (i.e. from Feb. 2004 to Nov. 
2010). The observations from the various sites are shown with 
different symbols. It is clear from Fig. [2] that during 2007 and 
2009 relatively dense observations over several months were 
carried out. Some variations in the amplitude of the brightness 
modulation can be seen over the six years. For example, dur- 
ing the first observing season » 0.5 mag fainter minima were ob- 
served while during the 2009/2010 observing season « 0.3 mag 
brighter phases were recorded. Fig.[3]shows the light-curve dur- 
ing the 2009/2010 observing season. The arrows show the ex- 
pected maximum positions based on the period measured (see 
Section 0J, which agree very well with the observed maximum 
positions. Variations in the peak brightness by up to « 0.5 mag 
are evident from Fig. [3] Very interesting is the enormous bright- 
ening of CHS 7797 at JD 2455214 with I » 16.7 mag, i.e. more 
than 1 mag brighter than the typical peak values. We will discuss 
the possible nature of this unusual "supermaximum" in Section 

m 



2004 and 2009 with the 2m Himalayan Chandra Telescope 
(HCT) at the Indian Astronomical Observatory and the 2.3m 
Vainu Bappu Teles cope (VBT) at the V ainu Bappu Observatory, 
both in India (see Pari har et alJ 12009. for details). The com- 
bination of all these data provides a unique I-band data set 
consisting of 663 observations over a total of seven observing 
seasons. Additional optical observations at other wavelengths 
were obtained at the MMO (R and z') and from HCT and 
VBT (R and V). In addition, J, H, K s images are available 
from observations with the CFHT Wide field InfraRed Camera 
(WIRCam) and the UKIRT Wide Field CAMera (WFCAM) 
carried out in 2006 and 2009, respectively (Alves de Oliveira, 
priv. comm.). We complemented our data set with recently re- 
leased Spitzer/IRAC observation of CHS 7797 using the IRAC1 
(3.6jum) and IRAC2 (4.5//m) filters taken as part of the YSOVAR 
project during 2009 (Morales -Calderon et al .1201 lh . Table 1 pro- 
vides an overview of the multiwavelength photometric observa- 
tions on CHS 7797. Very important for a detailed monitoring of 
CHS 7797 were the Spitzer observations, the WFCam observa- 
tions, the WIRCam/CFHT observations, and the I-band obser- 
vations (MONET/N, MMO, and HCT) between October 23 and 



4. Period determination 

The period determination is based on the I-band observations, 
since they provide the by far most complete and dense data set, 
extending over a very long time base. In addition, young very 
low mass objects are brighter in this part of the spectrum and, 
particularly important, I-band images suffer much less from neb- 
ular background contamination than R-band images. 

To determine the periodicity of the light modulation we 
used a combination of two pe riodogram ana lysis techniques: 
the Lomb-Scargle periodogram (Scargle1 [T982T) and the CLEAN 
algorithm dRoberts et al. 1 1987t> . As described in detail in 
Rodriguez-Ledesma et al (2009), these two techniques deliver 
a reliable period determination in which the effects of alias- 
ing and false peaks due to inhomogeneous sampling no longer 
appear. We derived false-alarm probabilities (FAP) based on 
the Lomb-Scargle periodogram to test the significance of the 
period found. Independent of the Lomb-Scargle periodogram 
FAP calculation, we used the statistical F-test and a de rived 
FAP from it, as described in | Scholz & Eisloffell (l2004bl) and 
iRodrfguez-Ledesma et al.l d2009l) . Since the FAPf->«.»> represents 
the probability that the period found is caused by variations in 
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Fig. 3. 1-band light-curves for the 2009/2010 observing season during which observations were performed almost continuously over 
«: 4 months. Arrows show the observed and expected positions of the maximum light. During these season a "supermaximum" was 
observed at JD = 2455214 (see section[6]for details). 



the photometric noise, it is completely independent of the peri- 
odogram analysis. 

We determined a period for the individual seven observing 
seasons and for the whole data set from 2004 to 2010. The results 
agree within the errors in all cases. The derived period based on 
all 663 I-band observations is P — 17.786 + 0.03 days. The 
computed FAPs, based on the Scargle periodogram and the F- 
test, give extremely low probabilities (<1 x 10~ 15 ) that the de- 
tected period is not real. Owing to the very long time base of 
our observations, we carefully searched in the periodograms for 
additional lower power peaks corresponding to longer periods, 
but no other significant peak was found in either the Scargle or 
the CLEAN periodogram. 

Fig. [4] shows in the top panel the phased light-curve based 
on the whole I-band data set. The bottom panel of Fig. [4] shows 
the phased light-curves of the individual seven observing sea- 
sons. Slight differences in the shape of the phased-light curves 
for different seasons are clearly seen. For example, the phased 
light-curves for the years 2004 and 2005 show pronounced max- 
imum and minimum brightness phases (i.e. triangular shape), 
while during the seasons 2006/2007 and 2008/2009 the maxima 
of the phased light-curves are slightly flatter. It is also evident 
that there are some changes in the amplitude between seasons. 
A more detailed description and discussion of these features is 
given in Section|6] 



5. Multiwavelength data 

In addition to the optical data set, we analysed NIR observa- 
tions taken during November 2009 in the J, K s , IRAC1 (3.6/vm) 
and IRAC2 (4.5/im) bands. Fortuitously, some of the K s -band 
observations from WIRCam/CFHT, those in the J-band from 
WFCAM/UKIRT and I-band observations were taken simulta- 
neously with the Spitzer observations. These particular obser- 
vations are extremely valuable because they allow us to study 
colour changes of CHS 7797 as a function of wavelength and 
the spectral energy distribution (SED). In FigfSJ these SED^] in- 
dicate that CHS 7797 is a very red and embedded pre-main se- 
quence star. If we take into acco unt the classification based on 
the IR slope bv lLada et all J2005) . with a slope of 0.8, CHS 7797 
classifies as a Class I object. As evident from Fig. [5] the bright- 
ness difference between maximum and minimum light is nearly 
constant until 2.2 //m, and for the longer wavelengths the differ- 
ence is smaller. 

The time interval over which we have obtained simultane- 
ous observations in I, J, K s , 3.6 and 4.5^/m is displayed in Fig. [6] 
This time interval of ss 40 d covers more than two full period cy- 
cles and is best sampled by the IRAC data. The ptp amplitude 



1 We note that we used the simultaneous data for the I, J, K„ 
IRACl[3.6yum], and IRAC2[4.6 y um] bands along with the R, z', H, 
IRAC3[5.8;um], and IRAC4[8.0yum] data, which were taken at different 
epochs. 
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Fig. 4. Top: Light-curve of the whole I-band data set from 2004 
to 2010 phased with a period of 17.8 d. Bottom: Phased light- 
curves of the individual seven observing seasons shown with 
different symbols. Slight differences in the shape of these phased 
light-curves are visible. 
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Fig. 5. Spectral energy distribution (SED) of CHS 7797 dur- 
ing maximum (circles) and minimum (open squares) brightness. 
Open triangles show IRAC measurements at 5.6 pm and 8 pm, 
which are single-epoch observations taken about 3-4 days after 
maximum. The best-fitted model found by VOSA, i.e. NextGen 
models with T e // = 2700 K and log g — 4, is shown. Note 
that brightness difference between minimum and maximum de- 
creases for A > 2 /urn. 
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Fig. 6. Light-curves in I, J, K s , 3.6/mi and 4.5yum bands during 
part of the 2009/2010 observing season. 



0.6 in 4-.5/J.. Typical amplitudes in the observed bands are listed 
in Table|2] Furthermore, the shape of the light-curves between 
the I and IRAC bands clearly looks different. While in the for- 
mer case the light-curve is more sinusoidal with a pronounced 
minimum, the light-curves at the latter two wavelengths look 
more like brightness maxima superimposed on a roughly con- 
stant light level, resulting in practically flat and broad minima. 



show a slight decrease towards longer wavelength. The relative 
decrease is lower (20%) in the optical/NIR from 1.7 in I to 1.4 
in K.E than in the IRAC bands (35%), from 0.9 in 3.6jU and to 



2 Note from Table[2]ihat the ptp amplitudes in I and z' are the same 
and equal to 1.7 mag. 



5.1. Colour changes 

To study the colour variations as a function of wavelength in 
more detail, we display several colour and magnitude light- 
curves in Fig. [7] Although variability is clear for all photomet- 
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ric bands, not all colours displayed in Fig.[7]show variations that 
correlate with the magnitude variations, i.e. the object is not al- 
ways becoming redder or bluer at bright or faint states. 

To gain a better understanding of the colour-brightness re- 
lations in Fig. |7] we display the following three correlations 
in Fig. [8] (R-I) versus I, (I-z') versus I, and (J-K,) versus K. 
Moreover, to study such correlation at longer wavelengths we 
show in Fig.|9]the (K s -3.6/im) and (K s -4.5^m) colours as a func- 
tion of the K, band magnitude. 

As evident from Fig. [8] there is no correlation between 
colours and brightness (correlation coefficient R«0.2 in all 
cases), quite in contrast to Fig. [9] where a clear and sta- 
tistically significant correlation is observed between the K s - 
3.6//m and K s -4.5yum colours as a function of K, m agni- 
tude (R = 0.97). As addressed by ICarpenter etail d2001l) and 
lAlves de Oliveira & Casalil (120081) . the slope of the colour-flux 
correlation can be used to infer the origin of the variability. The 
colour-brightness correlations shown in Fig.[9]have in both cases 
a slope that agrees w ithin the errors with the extinction law from 
Cardelli et al.l(ll989l) . The e xtinction values at 3.6 ^/m and 4.5//m 
for R = 3.1 were taken from lChapman et ail d2009l) . These clear 
correlations indicate that the variability is probably caused by 
variations in the extinction. The lack of any colour magnitude 
correlation for A < 2 fim indicates that the absorbing dust parti- 
cles have typical sizes of 1-2// m and are therefore larger than 
the typical dust particles in the interstellar medium. Fig. [8] only 
suggests a weak tendency for CHS 7797 to becomes bluer when 
fainter for A < 2 yum, which if real, might be due to scattered light. 



Table 2. Relevant properties of CHS 7797. 



Parameter 


Value 


Position(2000) 


RA=05 35 08.7 DEC= -05 31 27.3 


Photometry " 


Maximum (JD=2455 137) 


Minimum (JD=2455 130) 


I 


17.35 


19.01 


R-I* 


2.25 


1.9 


I-J 


2.24 


2.25 


J-K 


1.4 


1.5 


K,-3.6 


1.3 


1.8 


K s -4.5 


2.2 


3.1 


Typical ptp amplitudes ' 


1=1.7 z'=1.7 J=1.6 K= 


1.4 [3.6]=0.9 [4.5]=0.6 


Spectral type 


M6<' 


T eff 


2840-3000 


Distance (pc) 


450'' 


Luminosity (L Q ) 


0.01-0.06' 


Mass (M Q ) 


0.055-0.1 1 


Age (Myr) 


1 





" when available, the quasi-simultaneous photometric observa- 
tions were used for the typical photometric values, since as seen in 
Fig.|3]there are variations in the minimum/maximum levels. 
* the JD maximum is 2455474 and minimum 2455468 in this case. 
' typical R-band ptp amplitudes are not listed because CHS 7797 could 
not be reliably measured at minimum brightness. 
d from Rodriguez-Ledesma et al 2012, in prep. 
' from lO'deli fcOOlh 
f from this work. 



5.2. Stellar and disc properties of CHS 7797 

We followed several approaches for estimating the luminosity 
and mass of CHS 7797. We used the fluxes at maximum light 
corresponding to JD 2455137 (see the maximum peak in Fig. [6] 
for which I, J, K„ IRAC1 and IRAC2 observations are avail- 
able). We also considered that CHS 7797 might be a binary sys- 
tem composed of identical spectral type (SpT) components (see 
Section[6]). We point out that the classical approach based on 
the dereddened J flux has to be considered with extreme cau- 
tion since we are not dealing with normal interstellar dust, but 
with larger dust particles. Therefore we focused our luminosity 
and mass estimates on methods that are not, or only slightly, de- 
pendent on extinction. 

From a preliminary spectral analysis, and based on four 
narrow-band spectral indices particularly selected to measure 
the change of strength of molecular absorption feautures with 
spectral type, we derive a SpT of M6+ 0.5 for CHS 7797 (for 
details see Rodriguez-Ledesma et al 2012, in prep.). W e used 
the SpT-T e yy relations derived by iLuhman et al.l (120031) . which 
is best suited for young objects. For an M6 young dwarf, the 
T e ff based on this relation is 2990 K. If the commonly used SpT- 
T e ff relation for fiel d dwarfs is use d, one derives for CHS 7797 
a T„ff of 2840 K dLuhmanl [l999). Based on the models of 
Bara ffe et al.l (1 19981) at 1 Myr, the corresponding luminosities 
and masses for T eff = 2840/2990 K are L = 0.025/0.065 L and 
M = 0.055/0.1 M Q . It is important to note that the SpT-T c// re- 
lations for very young objects are relatively uncertain, and at 
these ages the available models are still relatively inaccurate. 
Nevertheless, CHS 7797 probably has a mass and luminosity 
within the previously mentioned values. If CHS 7797 is a binary 
system with identical spectral type components, we expect that 
each of the stars has a mass between 0.055 and 0. 1 M . 

As a second approach , we used the available VOSA SED 
analyser (Bav o et alJ2008h to compute the best set of parameters 
from the SED at maximum light. VOSA performs a statistical 
test to decide which set of synthetic photometry best reproduces 
the data. During the fitting process, the tool detects possible in- 
frared excesses and does not consider those data points of the 
SED where the excess has been found, i.e. in our case the IRAC 
bands are not taken into account in the fitting. The best-fit corre- 
sponds to Teff=2700 log g=4 in the NextGen models, with aj 2 
of 2. The derived bolometric luminosity is L = 0.024 L Q if A,,=0 
is considered and L = 0.040L G if we assume an A,,=2.45 This 
luminosity is the lumino s ity of the system. Assuming a binary 
system, the iBaraffe et al.l (Il998l) models predict a mass of each 
component of 0.04 - 0.055 M , in agreement with the derived 
masses and luminosities based on the SpT-T e // scale. 

In addition to VOS A SED analyser, we m ade use of the SED 
fitting tool described in lRobitaille et al.l d2007l) . This SED fitting 
tool uses the 200,000 model SED s for young stellar objects pre- 
sented in Robitaille et al.l (120061) . The models are characterised 
by 16 stellar, envelope, and disc parameters such as stellar mass 
and temperature, the envelope accretion rate, disc mass, accre- 
tion rate, and inner disc radius. To infer some physical proper- 
ties of CHS 7797 and its disc, we compared the modelled SEDs 
with the observed SED taking into account the ten best-fitted 
models. The derived temperature of the central source ranges 
from 2700 to 3000 K and the mass between 0. 1 and 0. 1 8 M , for 
the different best-fitted models. It is important to note that the 

3 T his value is derived by using intrinsic (I-J) colours from Leggett 
( 1992) at maximum brightness. Nevertheless, we emphasise that the 
derivation of an A v has to be taken with caution since we are dealing 
with grey dust at A < 2 //m. 
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Fig. 7. Colour and magnitude variations for two different time intervals. The top-left panel shows the variations in R-I, I, and R. The 
other three panels show colours and magnitudes for wavelengths up to 4.5 fim. No colour variations are evident for the R-I and J-K 
colours, while strong variations are visible for A > 2 /an. In these later cases the star becomes redder when becoming fainter (see 
also Fig.[9j». 



minimum mass considered in the models is 0.1 M Q . All models 
suggest that the inclination angle of the system is between 75 
and 82 degrees, and therefore we are seeing the system almost 
edge-on. As stated above, in addition to the stellar parameters 
these models allow for the characterisation of the circumstellar 
disc. The derived disc mass for CHS 7797 ranges from 1x10 4 
to 5x1 0" 3 M Q . 

As shown in Fig. [5] and also as computed with the VOSA 
SED analyser, the NIR excess starts at about 3/an, which sug- 
gests the presence of a gap in the inner disc. Based on the SED 
fitting tool, the inner disc edge in CHS 7797 is most probably 
located between 0.2 and 0.6 AU. The total range is larger (min- 
imum of 0.04 and maximum at 2 AU) but most of the derived 
values are between 0.2 and 0.6 AU. It is interesting to compare 
the derived inner disc edge with the radius at which the Keplerian 
period corresponds to 17.8 d (R pe riod)- Assuming stellar masses 
of 0.07 and 0.2 M , R period is 0.05 and 0.075 AU, respectively. 
These values are much lower than the typical inner disc edge 
values derived. Only one of the best-fitted models suggests that 
the inner disc edge is slightly smaller than R pe riod- All other nine 
models indicate a much larger inner disc edge (up to 50 times 
larger than Rperiod), and therefore suggest the existence of an in- 
ner disc gap. 



To summarise, we find that CHS 7797 is a very low mass 
object whether it is a single or binary star, with masses that 
are probably substellar or close to the substellar boundary. In 
addition, there are indications of an inner disc gap. 



6. Discussion 

6.1. The nature of the "internal clock" 

The variability observed in CHS 7797 is undoubtedly periodic 
with a period of 17.8 d. The periodic flux variations in young 
stars are usually interpreted by rotational modulation attributed 
to starspots in their atmospheres. The typical light-curves pro- 
duced by these spots show quasi-sinusoidal shapes. In addition 
to a rotational motion, orbital motions are intrinsically periodic 
and can of course give rise to a periodic flux variation. As out- 
lined in the introduction, this is true not only in eclipsing binary 
systems, but can also be caused by periodically variable absorp- 
tion in a circumstellar/binary disc. 

The 17.8 d period measured for CHS 7797 is rather long 
compared to the typica l rotational periods measured in young 
very low mass objects. Rodriguez- Ledesma et aD (|2009) found 
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ing that any variation due to variable extinction by dust particles 
similar to those found in the interstellar medium would cause 
a quite different colour magnitude correlation. The dotted lines 
show the least-squares fit to the data, which is not statistically 
significant in any of the three correlations. 



a median rotational period of a; 2d for 139 objects with esti- 
mated masses < 0. 1 M in the ONC. In the same work no objects 
with rotational periods > 13 d were reported at these low masses. 
Studies in other young clusters reveal similar values for the rota- 
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ison. The positive slope in both correlations agrees with variabil- 
ity caused by variations in extinction by dust particles smaller 
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tional periods (e.g. Joerg ens et alJ2 003, found rotational periods 
of 2-3 d for brown dwarf candidates in Cha I). 

Since co ol spots can only p roduce brightness modulations 
I < 0.7 mag dHerbst et alj 1 1994b . the very large amplitudes of 
la; 1.7 observed in CHS 7797 can only be caused by rotational 
modulation if we are dealing with a hot spot. If the light mod- 
ulation of CHS 7797 were caused by a hot spot, Fig. [8] should 
show a correlation between brightness and colour, i.e. the star 
should become redder when becoming fainter. This is not ob- 
served. In addition, the variation amplitudes caused by hot 
spots must be larger in the optical than in the NIR, particu- 
larly for a cool object like CHS 7797, because the relative con- 
tribution from the cool stellar atmosphere is strongly increas- 
ing towards the NIR, while the relative contribution from the 
hot spot is correspondingly decreasing. This is not observed 
for CHS 7797, since the amplitudes do not decrease with wave- 
length for A < 1 yum, a nd only slightly for A ^ 1 /im. Following 
the same approach as ICarpenter et alj (1200 ll) . we investigated 
the amplitudes of the brightness modulation that can be theo- 
retically produced by hot and cool spots. We fixed the stellar 
T e ff to 3000K, and assumed that spots can be characterised by 
single-temperature blackbodies covering a certain fraction / of 
the photosphere. We investigated the amplitudes at R, I, J, K s , 
3.6pm, and 4.5//m, for hot spot temperatures of 5000 K, 6000 K, 
7000 K, 8000 K, 9000 K, and 10000K and for hot spots cov- 
ering 1%, 5%, 10%, 20%, and 30% of the stellar surface. We 
find that to produce the amplitudes observed in CHS 7797 at 
A > 1 fim (AJ = 1.6 ,AK. S = 1.4 ,A3.6/mi = 0.9 ,A4.5/mi = 0.6), we 
need very hot spots (TalOOOOK) and spotted areas covering 
more than 20% of the stellar photosphere. These bright hot spots 
covering such a large part of the stellar surface would produce 
amplitudes in R and I bands of a 4 and w 3 mag, respectively, 
which are much larger than the a; 1 .7 mag observed in CHS 7797. 
In addition, variability studies of TTS with hot spots often show 
phase jumps and substantial change s in their modulati on am- 
plitudes from epoch to epoch (e.g. iHerbst et al.l 1 19941) . while 
we observed only a moderate change in amplitude over the six 
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years. We performed the same test for cool spot temperatures of 
1500 K, 1700K, 1900K, 2100K, 2300 K, 2500 K, and 2700 K 
(i.e. 10% to 50% cooler than the photosphere) and for cool spots 
covering fractions from 10% to 75% of the stellar surface. We 
found that not even the coolest 1500 K spot covering 75% of the 
surface of CHS 7797 can produce the large amplitudes observed 
in R, I, z, J, and K bands. 

We conclude from the previous arguments that a rotational 
light modulation due to hot or cool spots cannot explain the pe- 
riodic variability of CHS 7797, and therefore we most probably 
are dealing with an orbital motion. 

6.2. Nature of the variability 

The most likely source of orbital variability in an object such as 
CHS 7797 is variable extinction from circumstellar matter, most 
probably located in a disc. 

Our multiwavelength observations demonstrate that 
CHS 7797 shows no correlation betw een colour and fl ux for 
A<2ym. The NIR monitoring by ICarpenter et alj d2001l) 
also clearly indicates no colour-flux correlation in the NIR 
J,H,K;s bands for CHS 7797. Interestingly, we found a strong 
colour-flux correlation at A > 2 /mi that agrees with variability 
caused by variable extinction (see Fig.|9j>. The lack of a colour- 
magnitude correlation at A<2fim suggests that grain-growth 
has happened among the occulting material. Grain sizes with 
a minimum size of about 1 - 2yum in the circumstellar envi- 
ronment can account for the grey absorption at A < 2 fim, while 
the distribution of grain sizes > 2 /im should be similar to the 
interstellar distribution to explain the colour-flux correlation 
found at longer wavelengths, as clearly seen in Fig. [9] We note 
that various studies in the last decades s howed that discs in TTS 
usually undergo grain processing (e.g iPrzygodda et alJ [2003; 
iKessler-Silacci et al.1 120051: IShuping et alj 120061) . Based on 8 
- 1 3 low-r e solut ion spectra of eight proplyds in the ONC, 
Shuping et al. (2006) found that the typical grain size in the 
surface of the discs is a few microns, similar to what was found 
in young TTS in other star-forming regions. 

6.3. Is the absorbing dust located in a circumstellar or 
circumbinary disc? 

CHS 7797 could be a single object with a circumstellar disc and 
an structure or protoplanet orbiting with a period of 17.8 d or it 
could be a low-mass stellar or brown dwarf binary system similar 
to KH 15D. 



6.3.1. Single-object scenario 

In the single-object scenario, the eclipses in CHS 7797 are 
caused by variable amounts of circumstellar dust in the line-of- 
sight to the object. Since CHS 7797 is at the maximum bright- 
ness state only for « 1/3 of the orbital phase, the structure in 
the circumstellar disc causing the ocultation has to be quite ex- 
tended. This structure in the circumstellar disc is expected to 
orbit CHS 7797 with a period of 17.8 d. This scenario can only 
be taken into account if there is no inner disc cavity comparable 
in size with R period- A large inner disc gap, as suggested by nine 
of the ten best-fitted SED models, would make the single model 
scenario rather unlikely. Nevertheless, since one of the ten best- 
fitted SED models indicates that the inner disc edge is < R per iod, 
the single object scenario cannot be completely excluded. 



In the following, on the basis of the observ ations of V7 1 8 Per 
dGrinin et alJ200a 12008: Her bst et al.l2010bh . we argue that sin- 
gle periodically variable PMS stars with large variation ampli- 
tudes can be explained by strong absorption in their circumstel- 
lar discs over large fractions of their orbit. 

The TTS V718 Per is located in the 3Myr-old cluster IC 348. 
It shows periodic variability by 0.8 mag in I with a period of 
4.7 yr and eclipses that last « 3/4 of the orbit dGrinin et al.l l2006. 
2008). For comparison, the eclipses in CHS 7797 l ast 2/3 of the 
orbit. The light-curve of V718Per (see Fig. 1 in iGrinin et al.l 
2008) looks smoother than the light-curve of CHS 7797, but 
they are similar in shape. Although in the first attempts to ex- 
plain its unusual periodic variability a binary system similar 
to KH 15D was consid ered (see also next section), radial ve- 
locity measurements by IGrinin et al.l (120081) exclude the binary 
model for V7 1 8 Per and put an upper limit to any companion to 
~6Mj up . The authors concluded that the star is surrounded by 
a nearly edge-on circumstellar disc with an irregular mass dis- 
tribution presumably at the inner disc edge. As pointed out in 
their work, the extended circumstellar feature could be an irreg- 
ular azimuthal mass distribution in the disc, such as spiral arms 
due to the interaction of the disc with a pla net/protoplanet, or a 
warped disc structure. iHerbst et al.l d201 0b) considered a similar 
model in which the disc structure is strongly altered by a massive 
planet/protoplanet. 

An additional exam ple of this kind of variability is presented 
in IScholz et ail (120091) . Based on a NIR monitoring of 30 low 
mass objects in crOri (3Myr), IScholz et al.l (|2009) found evi- 
dence of a very low mass object (i.e. star #2 in their work) in 
which variability might be caused by inhomogeneities at the in- 
ner edge of the disc. In all cases, the system has to be seen almost 
edge-on to account for the extinction variations. 

AA Tau-like stars also show variability due to inner disc per- 
turbations that are associate d with a magnetic interaction be- 
tween the star and the disc (Ale ncar et al. l2010l) . AA Tau-like 
light curves are found to be fairly common among young stars, 
but the comparison with the light curves of CHS 7797 reveals 
significant differences. First, AA Tau-like stars are only quasi- 
periodic, (i.e. show time evolution over a few cycles) and not at 
all as perfectly periodic as CHS 7797 (FAP = 10~ 15 ). Second, the 
light curves of AA Tau systems have a flat maximum that is in- 
terrupted by minima events that are highly variable both in width 
and depth, quite different from the shape of the light curves of 
CHS 7797, as shown in Fig|3]and|4] 

6.3.2. Interpreting our data in terms of a KH 1 5D type binary 
model 

In the binary scenario, CHS 7797 is a low-mass stellar or brown 
dwarf binary system with an orbital period of 17.8 d and a 
circumbinary disc. The orbital period of 17.8 d and the esti- 
mated separation between the binary components of ^ 0.08 AU 
(for Mi = M2 =0.1 M ) are consistent with other known young 
low-mass binaries (e.g. St assun et al.l 12006). By analogy with 
KH 15D, we assumed that the disc is slightly inclined with re- 
spect to the orbital plane of the binary and that only one of 
the components can be seen during apastron through a more 
transparent outer region of the disc during a 1/3 of the orbit. 
This means that - in contrast to KH 15D, for which until 2007 
one star (star A) was unocculted by the circumbinary disc dur- 
ing maximum brightness - CHS 7797 is always occulted by the 
disc but the amount of extinction is periodically changing de- 
pending on the position of one of the stars with respect to the 
more transparent outer regions of the circumbinary disc. Note 
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that the « 1.7 peak-to-peak amplitude of CHS 7797 in the I band 
is much smaller than the nearly 4 mag amplitude observed in 
KH 15D. One possibility to account for the variations at maxi- 
mum level in CHS 7797 and the general variability in I-band of 
typically 0.3 mag at a given phase, is a disc that does not ex- 
tend smoothly and homogeneously to its outer edge but that is 
highly structured, and therefore we see the star at the same or- 
bital phase through portions of the disc that are correspondingly 
variable in transparency (in some cases from one orbital period 
to the following one). These variations in transparency can in 
principle also explain the "supermaxima" that we observed in 
our data (see Fig. [3}, in which the object is about 1 mag brighter 
in I band than the typical maxima level. Over the six years of 
monitoring, only two supermaxima were detected (an additional 
one with large photometric errors was ignored), which clearly 
indicates that these extreme increases in brightness are rare. We 
note that the "supermaxima" events do not occur at maximum 
phase, but some days after maximum (« at the end of ingress 
phase). Another possible explanation for the observed "super- 
maxima" events is flaring activity, which cannot be excluded for 
CHS 7797 because it is a strong Ha emitter (Rodriguez-Ledesma 
et al., in prep.). For the same reason the » 0.3 mag variations 
can be accretion-related. We expect the latter effect to contribute 
only in a; 0.05 mag colour changes in the (R-I) versus I dia- 
gram (Fig. [8]i, and we note that some CTTS s show little colour - 
brightness correlations (e.g. AA Tau, DR Tau lHerbst et al. 1994). 

The present interpretation of the KH 15D system considers 
the precession of the circumbinary disc, which causes the occult- 
ing sharp edge to move ac ross the orbit of the binary dWinn et al.l 
2006: Herbstetal.ll2 010a). This substantial difference can be ap- 
preciated if one compares the light-curves of CHS 7 797 in Fig. [3] 
with t he light-curves of KH 15D (see e.g. Fig. 3 in iHerbst et alj 
2010a): KH 15D light-curves have a constant maximum bright- 
ness within the same observing season, and only show a steady 
and smooth decrease in amplitude with years. Furthermore, the 
shape of the light-curves in KH 15D also changes in a continu- 
ous and smooth fas hion from one obse rving season to another. 
As seen in Fig. 3 o flHerbst et al](l2010al) . the star is at bright state 
for a shorter amount of time every consecutive observing season. 
Over the six years of monitoring of CHS 7797 we did not ob- 
serve such a monotonic decrease (or increase) of the maximum 
light duration. This implies that any possible disc precession in 
CHS 7797 occurs at a much lower rate than in KH 15D. We note 
that the typical shape of the light-curve brightness modulation 
shown by CHS 7797 is similar to that shown by KH 15D during 
the 2007/2008 season (Fig. 3 in lHerbst et alll2010al) . 

7. Conclusions 

CHS 7797 is undoubtedly a very interesting and unusual young 
and possibly substellar object that is occulted by circum- 
stellar/circumbinary matter. Whether it is a single object or 
a binary system, CHS 7797 is surrounded by a circumstel- 
lar/circumbinary disc that is viewed almost edge-on and in which 
substantial grain growth has taken place. A comparison of the 
light-curves and phased light-curves of CHS 7797 with the three 
other known "unusual" periodic variables (i.e. KH 15D, WL4 
and V718 Per) suggests some similarities but also discrepancies. 
In particular, the other three objects show much smoother light 
curves than CHS 7797, which implies that, for as yet unknown 
reasons, the disc around CHS 7797 is much more inhomoge- 
neous, and/or that there is a contribution of accretion-related ir- 
regular variations. Clearly, CHS 7797 deserves further attention. 
An upcoming spectroscopic analysis (Rodriguez-Ledesma et al., 



in prep) will try to provide additional constraints on the nature of 
CHS 7797. Nevertheless, since CHS 7797 is very faint and since 
there is a strong background contribution from the Orion nebula, 
it will be hard to search for radial velocity variations to confirm 
the possible binary character of the system even with the instru- 
mentation currently available at 8-10m class telescopes. 
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Appendix A: Absolute photometry of stars in the 
CHS 7797 region 

We have performed Bessell R, I, and SDSS z' absolute photom- 
etry of nine stars in the vicinity of CHS 7797 based on obser- 
vations at MMO. These stars were used as secondary standards 
for the relative photometry described in Section 3 and have b een 
selected from the work of Rodriguez- Ledesma et al.l (2009) as 
non-variable (or very low variable) stars. Although the ONC is a 
well-studied region, due to the bright and inhomogeneous neb- 
ular background no absolute R and z' photometry was available 
so far. We therefore performed the absolute photometry by ob- 
serving Landolt standard fields and Sloan Sky Digital Survey 
(SSDS) fields. The table below summarises the results. 



Table A.l. MMO R, I, z' photometry for ONC stars in the vicin- 
ity of CHS 7797. 



ID" 


R 




I 




z' 


z ' 

*■ err 


255 


14.643 


0.071 


14.019 


0.071 


15.364 


0.084 


209 


15.748 


0.071 


14.315 


0.071 


14.727 


0.084 


10318 


15.557 


0.100 


14.285 


0.095 


14.798 


0.110 


148 


16.933 


0.089 


15.304 


0.082 


14.923 


0.105 


105 


18.258 


0.145 


16.015 


0.108 


15.144 


0.122 


120 


16.304 


0.071 


14.572 


0.070 


14.524 


0.084 


132 


16.342 


0.077 


14.706 


0.071 


14.734 


0.084 


10263 


18.577 


0.158 


17.028 


0.102 


17.269 


0.207 


103 


18.432 


0.138 


16.350 


0.088 


15.553 


0.130 



" ID numbers are from lHillenbrandl d!997h . 



